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Abstract 

Kinetics of silicon dry oxidation are investigated theoretically and experimentally at low tem- 
perature in the nanometer range where the limits of the Deal and Grove model becomes critical. 
Based on a fine control of the oxidation process conditions, experiments allow the investigation of 
the growth kinetics of nanometric oxide layer. The theoretical model is formulated using a reac- 
tion rate approach. In this framework, the oxide thickness is estimated with the evolution of the 
various species during the reaction. Standard oxidation models and the reaction rate approach are 
confronted with these experiments. The interest of the reaction rate approach to improve silicon 
oxidation modeling in the nanometer range is clearly demonstrated. 
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I. INTRODUCTION 



An important challenge imposed by CMOS downscaling is the growth of ultra-thin ox- 
ides of silicon of high quality, with a tight thickness control and a good uniformity jl]. Very 
important progress have been done during the last decade to improve the structural prop- 
erties and electrical limits of ultra-thin (< 4 nm) silicon oxide and oxynitride [2j. On the 
other side, the modeling of ultra-thin silicon oxide growth remains a difficult issue for the 
microelectronics industry. 

To date, the seminal work of Deal and Grove remains the main approach used in process 
simulators to describe the oxide growth [3]. Assuming a steady state reaction between 
molecular oxygen and silicon, Deal and Grove deduced that the oxide growth can simply be 
described by the following equation : 
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A 2 /AB ' (1) 
where the oxide thickness X is described by a linear-parabolic relationship as a function of 
the oxidation time t. The term ■§ characterizes the initial linear rate growth. The parabolic 
rate constant B governs the diffusion limited regime. An initial time offset r is necessary 
to take into account the presence of an initial thick oxide layer. The success of this model 
can be explained by its ability to describe both dry and wet oxide in the case of thick oxide. 
Another major achievement is the introduction by Kao et al. 4j of the influence of strain 
effects in the Deal and Grove approach. Using this generalized Deal and Grove approach, 
Kao et al. explained the retardation effect observed for the oxidation of curved surfaces jlj] . 

However, this model suffers from strong limitations such as the failure to describe the fast 
initia. growth regime for dry oxidation Q or the d.fflcu.ties to extend it to oxynitridation 
process [7|]. Much research work has been dedicated to elucidate the breakdown of the Deal 
& Grove model in the ultra-thin regime. The idea was to correct the Deal and Grove model 
by adding some empirical terms to the growth rate 8|, |9j. However, the physical origin of 
these terms is subject to discussion. 

An alternative approach has been proposed by Wolters et al. who considered that mo- 
bile ionic species are responsible of the silicon oxidation process 10|, 111]. Ionic transport 
and induced electric field considered in this model could explain some of the experimental 
phenomena observed in silicon oxidation such as the rate dependence on the orientation or 



the cross-over effects 12|. But the role of charged species in the oxidation mechanism is still 
debated even, for more fundamental simulations such as density functional theory [13I ll^j ]. 
Recently, a new alternative linear parabolic rate equation has also been proposed by Watan- 



abe et al. 



15| without the rate-limiting step of the oxidation reaction but assuming that 



diffusivity is suppressed in a strained oxide near the silicon interface. 

From the above discussion, a revision of the Deal and Grove model is clearly necessary 
to improve process simulation. An ideal model for the microelectronics industry must give 
access to : (i) a fundamental and quantitative understanding of the interaction between 
oxygen , silicon and the standard dopant at the Si/Si0 2 to better describe segregation effects 



16 



171 ] or to estimate the dose- loss of dopant, (ii) the amount of silicon interstitials at the 



interface (iii) the amount of nitrogen incorporated during a more complex oxynitridation 
process (NO, N 2 gas) to estimate the gate leakage current |18 |. The development of 
such model still represents a big challenge. Moreover, a multiscale modeling of silicon 
oxidation [jjj will probably be mandatory to simulate the growth of ultimate oxides ( ~ 1 
nm) fabricated by advanced oxidation process such as plasma nitridation step [h]]. 

In this paper, the modeling of dry oxidation kinetics at low temperature in the nanome- 
ter range is investigated. This paper is organized as follows. In section II, the reaction 
rate formalism is presented. The calibration method and some preliminary results are also 
discussed. Experiment at low temperature are described in section III. The two standard 

n nn 

models (Deal and Grove [3|], Massoud [6|, |8[) and the reaction rate approach are confronted 
with the experimental data in section IV. Most important conclusions are summarized in 
section V. 



II. THE REACTION RATE APPROACH 



Introduction 



An interesting study in silicon oxidation modeling has been published by Almeida et 

al. 
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21] . The main idea of this work is to analyze the main assumptions of the Deal 
and Grove model and to propose a more rigorous approach. The first approximation made 
by Deal and Grove is to consider that the reaction strictly takes place at the interface. 
The second approximation is the steady state regime which imposes the a balance between 
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oxidizing species entering the SiC>2 surface, the number of molecules diffusing into the oxide 
and the molecules reacting at the Si/Si02 interface. Thus the concentration of the oxidant 
species is not time dependent. As outlined in an very elegant way by Almeida et al., these 
assumptions are certainly the key to explain the limitations of the Deal and Grove model in 
the ultra-thin regime for dry oxidation. Moreover these assumptions limits the development 
of more sophisticated treatment of the oxidation mechanism. 

B. The model 



In this paper, a reaction rate approach similar to that of Almeida et al. [20|, |2l| has 
been adopted for the modeling of our experimental kinetic at low temperature. The main 
diffusing species is assumed to be molecular oxygen (O2) as demonstrated by several isotopic 
experiments [22j. The two other species are pure silicon (Si) and silicon dioxide (SiO^). The 
reaction of oxidation is simply given by: 

K 

Si+ 2 Si0 2 

where K is the reaction rate. Keeping in mind that the investigation focuses on the nano- 
metric regime, a [100] silicon surface associated with a one dimensional system of coordinates 
is considered (see Figure [1]) . The area to be oxidized is viewed as an assembly of silicon 
monolayers. The grid is defined to match as fine as possible this assembly of planes. This 
means that the vertical mesh step corresponds to the distance between two silicon planes 
(1.35 A). The oxygen flux in the furnace is perpendicular to the silicon interface. Since the 
model involves an oxidation mechanism in the nanometer range, the growth of the film is 
restricted to the vertical direction. 

The relative concentration rij of the species j (where j— Si, 2 , and SiO^) is defined as 
follow : 



n,(x,t) = ^j^ (2) 



where CWt) corresponds to the concentration in the plane in units of number of atoms 
per unit surface, and C® is the maximum possible concentration in the plane. rij(x,t) can 
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be viewed as a layer coverage. For example, a value of 1 for nsi, means that the coverage of 
the layer is complete and corresponds to a concentration of 0.91 xlO 15 at/cm 2 of silicon. 

In order to setup the system of equations, the evolution of each species is described. 
During the oxidation process, the molecular oxygen diffuses in the silicon dioxide, reaches 
the silicon and reacts with it. Considering the diffusivity of the molecular oxygen D and 
the reaction rate constant K, the evolution of molecular oxygen can be given by: 

^ = V(D ■ Vn 02 ) -K-no 2 - n Si (3) 

In this expression, the reaction region is defined by the overlap between the oxygen and 
silicon species. Next, silicon consumption has is described by: 

—— = -K ■ n 02 ■ n Si (4) 

To sim plify the previous equation, the fact that some silicon interstitials generated during 
oxidation |23| has been neglected. The conservation law of the different species allows the 
description of the silicon dioxide creation: 

dt = K -n 0% - n Sl (5) 

The mathematical formulation of the three steps (diffusion of molecular oxygen, reac- 
tion with silicon, and creation of silicon dioxide) leads finally to a system of three coupled 
equations (Eqs [31 HI [5]). This system is numerically solved using a Cranck-Nicolson scheme 



24]. 



C. Boundary conditions 

To complete the mathematical description of the oxidation mechanism, boundary condi- 
tions need to be refined. For oxygen, it is necessary to evaluate the concentration of oxygen 
molecules that lies on the silicon surface. In this study, this concentration is estimated by 
simple physical considerations. 

Considering that the oxygen is described by the ideal gas law, the number of oxygen 
molecules in the furnace chamber can be expressed as a function of the gas pressure P and 
temperature T. Assuming that an ideal monolayer of oxygen molecules lies on top of the 
wafer surface, the density is : 
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N 02 NAPh 

C o 2 = — = -^r- (6) 



where Na is the Avogadro number, R the universal gas constant, S is the surface of the 
lattice and h the height of the oxygen layer. If the maximum density Cq^ x occurs when all 
the molecules form a square lattice of a = 3 A spacing: 



C™ x = - — ~ 1.0 x 10 15 at/cm 2 (7) 

*J mm 



The oxygen layer coverage on the surface can be expressed by the following expression : 



= = N^hP = P 

L>2 £imax RT T 



which further simplifies to the last expression of Eq. El if the height of the layer is about 
h = 3 A. 

The layer coverage of oxygen at the interface depends on both the partial pressure and 
the temperature. Figure [2] represents the variation of this coverage with the temperature. 



D. Calibration 

A rough estimation of the the oxygen concentration at the interface can be deduced using 
simple physical arguments. In order to estimate the value for the reaction rate K and for the 
diffusivity D, a calibration step has been undertaken on experimental kinetics. The large 
experimental oxidation database of Massoud (measured with an in-situ ellipsometer) has 



been chosen |6]. A numerical optimization method based on simulated annealing has been 
undertaken 25j. This method explores the whole solution space for the parameters (here 
K, D) to reach a global minimum for the error function. The probability of jump between 
different parameter values is proportional to Boltzmann probability distribution. At the 
beginning of the calibration step, the temperature Tg (for the Boltzmann distribution) is 
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high to explore all the solution space and then Tb is cooled down to limit the variation of 
the parameters in order to reach a minimum. The calibration procedure used is : 



• Start with an initial value of K and D and with an important value for the control 
parameter Tb 

• Calculate the theoretical oxide thickness X t h at the different experimental points 

• Estimate the error function E defined by the difference between the theoretical esti- 
mations and the different experimental thickness of the kinetics: 



Where Nexp is the number of the experimental points of the oxidation kinetics, a, are 
weighting coefficients defined such as to enforce the algorithm to minimize principally 
the error function for the lowest oxide thickness. 

• Modify the value of K and D to minimize the error function E. According to the 
Boltzmann probability associated with E, cool Tb- 

At the end of an important number of iteration steps (100000), the value of the two 
parameters minimizing the error function is obtained. 

E. Calibration results 

The parameters given by the optimization step to adjust the Massoud experimental data 
are reported in Table [B In order to extract a physical law for the variation of the two 
parameters, the different values have been calibrated again to match an Arrhenius law. As 
shown in figure [31 both the diffusivity and the reaction rate parameter agree well with an 
Arrhenius plot, which leads to the following law: 





D = D ■ exp — 



E D ' 



(10) 



k B T\ 
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The energy of activation is found to be Ed = 2.22 eV and the prefactor Dq = 1.291 x 10 11 
nm 2 /s. For the reaction-rate, we obtain similarly an Arrhenius law : 



K — Kq ■ exp 



E 



K 



k B T\ 



(11) 



with Kq = 2.022 x 10 7 s 1 and Ek = 1.42eV. It is worth noting that the activation energy 



26|. The 



and reaction rate are physically reasonable and close to published numbers 
calibrated activation energy for the diffusivity is in agreement with the experimental diffusion 
energy in silicon (2.42 eV) [26] and with ab-initio estimations (2.3-2.5 eV). The energy of 
activation (1.42 eV) for the reaction rate is not so different from the energy to break an Si-Si 
bonds (1.86 eV) (sT ] . 



F. Preliminary evaluation 

A primary consolidation of the present model has been performed through a comparison 
with published oxide thickness measured by ellipsometry. Figure H] compares the prediction 
of our modeling approach with the kinetics of Massoud et al. at atmospheric pressure from 
800°C to 1000°C. Each sample has an 1.0 nm initial oxide prior to oxidation. Thanks to the 
calibration step, an excellent agreement between the experimental data and the model is 
obtained. The oxidation rate for the different temperature is clearly well described both in 
the linear and the parabolic regime. No clear loss of predictivity is observed when using the 
analytical law for the reaction rate and the diffusivity rather than each pair of calibrated 
data. Next, the model is compared with the experiments made by Chao et al. 28J at 
atmospheric pressure in the nanometer range. The corresponding experimental data shown 
in the figure [5] were obtained with a multi-angle incident ellipsometer. Samples have a 
native oxide of 1.6-2nm covering the substrate. Similar conclusions can be drawn again 
considering that the variation of the oxidation growth rate is well described by the model 
for experimental conditions close to those used by Massoud et al. 

Finally, a comparison between our model and the original kinetics obtained by the Deal 
and Grove model has been performed. The main objective is to check the validity of the 
reaction rate approach to describe correctly the oxide growth in the thick regime. The 
comparison between the two models is of interest. Figure [6] presents the classical Deal and 
Grove kinetics simulated with the original parameters published in 1965 [3] . The breakdown 



8 



of Deal and Grove model to describe dry oxide growth below 20 nm is clearly emphasized. 
Almost no oxide growth is observed for a temperature near 800°C. The kinetics obtained 
by the reaction rate approach with the same temperature and pressure is represented in 
the figure [7J The differences between the two models is obvious in the ultra-thin regime. 
However the reaction rate approach is also able to describe the oxide growth in the thick 
regime. The growth rates predicted by the two models for various temperatures are in 
relatively good agreement. Even if the initial thickness considered in the two models is 
strongly different, the amount of oxide grown is almost identical a difference that does not 
exceed 5 % for 30 minutes at 1200°C. 



III. DRY OXIDATION AT LOW TEMPERATURE 



A complementary experimental work has been undertaken to validate the model. In order 
to reach the nanometer range with a conventional furnace, the oxidation is often performed 

nn 

by reducing the pressure of the growth ambient [29J, 130J , by diluting the oxidant gas with 



nitrogen 



3l| , or by lowering the temperature 32| . Oxidation with low thermal budget is an 



interesting solution. This point is illustrated by the study of Bhat et al. who propose the 
growth of ultra-thin oxides of silicon at low temperature [600°C-700°C] by wet oxidation 33 ]. 
In the present case, dry oxidation in a classical furnace is performed at low temperatures 
(725°C, 750°C) to fabricate oxide in the [1.5-4 nm] range. 



A. Oxidation experiments 

Low oxidation temperature is not sufficient to achieve oxide layers with interesting struc- 
tural properties in the nanometer range. Special care has to be given on the design of the 
thermal process to achieve a good homogeneity. Our experiments were performed on [100] 
silicon substrates (p type doped 5x 10 15 at/cm 3 ) that previously received a RCA cleaning 
and 1% HF dip. The oxide growth was carried out in a conventional furnace at atmospheric 
pressure. The oxidation process is divided into four stages as shown in figure M 

The sample introduction is performed at a temperature lower by 50° C than the main 
temperature step. Next, a temperature ramp is applied to reach the desired temperature 
during 30 minutes. The ambient is composed by a main flux of nitrogen (2 L/min) and a 
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small amount of oxygen (0.2 L/min). During this ramp, a very thin layer of silicon dioxide is 
created on top of the wafer before the start of the isothermal part of the oxidation step. An 
oxide of 1.53 nm (resp. 1.96 nm) is grown for a preoxidation ramp of 30 minutes at 725°C 
(resp. 750 °C) is created during this step. The thickness grown during the first part of the 
thermal cycle is nearly half of the maximum thickness range. However, this step is critical 
regarding the oxide thickness homogeneity. The preoxidation ramp has a clear impact on 
homogeneity which is less than 2 % on the whole wafer for the oxidation process at 725°C. 
Performing the same process with an inert ambient during the ramp-up leads to a strong 
increase of the oxide thickness dispersion up to 16 %. 

Finally, the oxidation stage is realized under an enhanced oxygen flux of (2 L/min). As 
shown in table II, the oxide growth is slow and the expected final thickness can be precisely 
controlled by the duration of the main oxidation stage since the reaction rate is very limited. 
The final step consists in an inert nitrogen atmosphere for a ramp down during 10 min. and 
by a thermal reflow of 30 min. No growth is expected during this phase. The main objective 
of the final process step is to improve the electrical properties of the oxide film and to reduces 
the amount of interface defects. 



B. Ellipsometry and TEM measurements 

Oxide thickness has been measured by a spectroscopic ellipsometer. In order to validate 
the ellipsometry measurements, some TEM analysis were performed. Two oxide layers 
grown with a preoxidation ramp at 725°C during 10 min and 30 min were analyzed. A 
polysilicon capping layer simulate the gate stack. The TEM analyzes (figure |8]) shows 
that a thin and uniform oxide layer. The following table (table [XT]) compare the thickness 
measurements given by the two methods. The good agreement between the ellipsometer and 
TEM measurements validate the ellipsometric method even if it has already been observed 



that for thickness below 2.5 nm, deviations become more important 



34)- 
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IV. SIMULATION OF OXIDE GROWTH AT LOW TEMPERATURE. 



A. Introduction 

The main objective of this section is to compare our experimental kinetics with the two 
standard models (Deal and Grove and Massoud) and the reaction rate approach. It has been 
previously verified that the reaction rate approach is able describe the kinetics in the case 
of thick oxides. However some questions still require clarification: i) is finally the reaction 
rate approach more adapted to describe the silicon oxidation in the nanometer regime ? ii) 
Is it possible to describe the oxidation in the low temperature regime where the oxidation 
reaction rate is very limited ? The objective of this section is to address these questions. 

On the other side, it must be kept in mind that oxide growth in this range [1.5-4 nm] 
is strongly dependent on processing parameters like e.g., the pre-oxidation ramp and sur- 



face cleaning methods [2lJ. Moreover, the experimental regime is relatively challenging for 
oxidation modeling. These two points emphasize the fact that only the oxidation kinetics 
variation during the main oxidation stage can be discussed. 

B. Comparison with Deal and Grove 

Our experimental kinetics data are compared with the two standard models. Since our 
objective is to tackle the model limitations, a continuous implementation of these models 
have been performed. The oxide thickness is obtained by integrating the following expression 
of the growth rate using very short timestep: 

d 4- = (12) 

dt 2X + A K 1 

In the present case, for the diffusivity and the reaction rate, calibrated parameters of 
standard TCAD tool has been used {35J]. This set of parameters is supposed to give more 
accurate results in the ultra-thin regime. Since, the oxidation temperature is below 1000°C, 
the linear reaction rate is govern by: 



— = 1.25 x 10 05 exp 
A 



1.76 eV 

k B T 



(nm/s) (13) 



and the diffusivity 
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B = 2.833 x 10 09 exp 



2.22 eV 
k B T 



(nm 2 /s) 



(14) 



Oxide growth has been simulated for the main oxidation steps at 725°C and 750°C. Since 
the initial oxidation ramp-up is not simulated, the theoretical curve must be shifted upward 
to match the experimental thickness obtained at the end of the temperature ramp up. As 
shown in figure [TOj the theoretical kinetic predicted by Deal and Grove [3J significantly 
deports from experimental data. The oxide growth rate is strongly underestimated. Almost 
no oxide is created at low temperature, in contradiction with our oxidation experiments. 



C. Comparison with the Massoud model 

As previously discussed, the Massoud model is an extension of the Deal and Grove model 
to improve the description in the ultra-thin regime. In practice, the growth rate expres- 
sion is corrected by adding a supplementary term that exponentially decays with the oxide 
thickness: C 2 exp( — j-) : 

dX B r ^ , X s , 

HT = 2xTa x[1 + c ^ p{ -T )] (15) 

C*2 involves with an Arrhenius law: 

2.56 eV 



C 2 = 61.52 exp 



k R T 



(16) 



and the oxide thickness of this additional term is controlled by the parameter length 
(L=7 nm). 

As shown in figure [TTj the introduction of the corrective term does not really improve 
the situation. Growth rate are still underestimated in the lowest temperature regime. It 
is worth noting that the situation improves at 750° C since it is much closer to the lowest 
temperature (800°C) used by Massoud et al. to perform their calibration. 
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D. Comparison with the reaction rate approach 

A final comparison has been performed between the experimental results and the reaction 
rate approach. Results are shown in figure O The first stage of the oxidation process is 
not considered since we are interested in simulating the main isothermal oxidation ramp. 
The theoretical oxidation kinetics are shifted upward in order to match the oxide thickness 
at the beginning of the oxidation ramp. Doing so, the agreement between the experimental 
kinetics and the model is fine specially for the oxidation at 725°C. the range of oxide thickness 
predicted by the reaction rate approach is clearly relevant. For 750°C, the model slightly 
overestimates the growth rate but the error remains reasonable (the maximum error is less 
than 0.4 nm for the complete set of experimental data). Probably the most interesting point 
is that the variation of the growth rate is relatively well defined. Moreover the oxidation 
kinetics seems to be not strictly linear as often described by the two standard models (the 
reaction rate approach predicts a behavior in oc t 8 ). 



V. CONCLUSION 



In conclusion, we have proposed a combined experimental and theoretical study for the 
dry oxidation of silicon at low temperature. A model based on the reaction rate approach 
and on recent developments [20] has been proposed. A set of parameters has been calibrated. 



A comparison with the two standard models shows a very good predictivity both in the thin 
and the thick regime. A complementary study has been carried out to test the predictivity 
of the two standard models and of the reaction rate approach in the nanometric regime 
([1.5-4 nm] range). Considering that only two parameters have been calibrated, the results 
obtained by this approach are very promising. Further work is needed to refine and improve 
the model. For example, it would be of interest to test if the reaction rate approach is able 



to describe the complex orientational effects observed experimentally by Irene et al. c 



or 



Ngau et al. [37]. A further extension of this approach could also be the modeling of various 
oxynitridation processes (NO, N2O) used by the microelectronics industry. 
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FIG. 1: Experimental system under investigation. Oxygen diffuses and reacts with silicon and 
creates silicon dioxide. A very fine ID mesh is defined to simulate the evolution of the various 
concentration. 
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FIG. 2: Variation of the layer coverage of oxygen with temperature. 
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FIG. 3: Reaction rate (K) and diffusivity coefficient (D) obtained by the calibration step for the 
different experimental kinetics of Massoud et al. for different temperatures. The Arrhenius laws 
are also reported. 
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FIG. 4: Experimental oxidation kinetics of Massoud (Points) et al. [6] and calculated with our 
model (lines). 
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FIG. 7: Oxidation kinetics predicted by the reaction rate approach for the same conditions (tem- 
perature, pressure). 




FIG. 8: XTEM analysis of an oxide layer of 19 A grown at 725°C during 10 minutes. 
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FIG. 9: The two thermal oxidation cycles at (725°C, 750°C) used in our oxidation experiments. 
Each thermal cycle is composed of a pre-oxidation ramp, an isothermal oxidation step and a thermal 
reflow. 
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FIG. 10: Experimental kinetics at low temperature for 725°C and 750°C (points) and the Deal 
and Grove model (lines). Oxide growth rate is underestimated leading to almost no oxide growth 
at 725°C. 
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FIG. 11: Experimental kinetics at low temperature for 725°C and 750°C (points) and the Massoud's 
model (lines). Oxide growth rate is clearly underestimated at 725°C. 
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FIG. 12: Experimental kinetics at low temperature for 725°C and 750°C (points) and our model 
derived by the reaction rate approach (lines). The first experimental kinetics is nicely described 
by the model whereas the second experimental kinetics is slightly overestimated. 



Tables 



T(°C) 




D (nm 2 /s) 


K (s- 1 ) 


800 


1.50e-04 


27.5 


7.0 


850 


1.43e-04 


124.4 


8.1 


900 


1.37e-04 


355.6 


15.8 


950 


1.38e-04 


1035.9 


23.8 


1000 


1.26e-04 


2268.1 


57.2 



TABLE I: The value of the diffusivity [D) and of the reaction rate (K) obtained by the optimization 
procedure to match the experimental kinetics of Massoud et al. [6J at atmospheric pressure (P=l 
atm). 



Duration 


Ellipsometry 


TEM 


10 Min. 
30 Min. 


18.90 ± 0.3 A 
27.10 ± 0.3 A 


19 ± 1.0 A 
27 ± 1.0 A 
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TABLE II: Validation of our ellipsometry measurements by TEM analysis on two samples oxidized 
at 725° C with the same pre-oxidation ambient N2/O2. 
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